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Abstract. We first give hypotheses of the bicharacteristic equations corresponding to the 
Enskog equation with an external force. Since the collision operator of the Enskog equation is 
more complicated than that of the Boltzmann equation, these hypotheses are more complicated 
than those given by Duan et al. for the Boltzmann equation. The hypotheses are very related 
to collision of particles of moderately or highly dense gases along the bicharacteristic curves and 
they can be used to make the estimation of the so-called gain and loss integrals of the Enskog 
integral equation. Then, by controlling these integrals, we show the existence and uniqueness of 
the global mild solution to the Enskog equation in an infinite vacuum for moderately or highly 
dense gases. Finally, we make some remarks on the locally Lipschitz assumption of the collision 
factors in the Enskog equation. 

P" 1 " 1. Introduction 



c3 



This paper is to consider the existence and uniqueness of the global mild solution to the 
Enskog equation with an external force in an infinite vacuum for moderately or highly dense 
gases. Throughout this paper, R + represents the positive side of the real axis including its 
origin and R 3 denotes a three-dimensional Euclidean space. In the presence of external forces 
depending on the time and space variables, the Enskog equation is as follows (see [7] or |10j): 

where / = f(t, x, v) is a one-particle distribution function that depends on the time t G R+, the 
position x G R 3 and the velocity v G R 3 , E(t,x) is a vector- valued function which belongs to 
R 3 and represents an external force with respect to the time and space variables, and Q is the 
Enskog collision operator whose form will be explained below. 

The collision operator Q is expressed by the difference between the gain and loss terms 
respectively defined by 
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d ■ Jn 3 xs[ 



Q+(f)(t,x,v) = F + (f)f(t,x,v')f(t,x-au;,w')B(v-w,uj)dujdw, (1.2) 



Q (f)(t,x,v)= / F (f)f(t,x,v)f(t,x + au;,w)B(v — w,uj)dojdw. (1-3) 
Jn 3 xs\ 

Here and below everywhere, = {to G S 2 : u>(v — w) > 0} is a subset of a unit sphere surface S 2 
in R 3 , a is a positive constant that represents a diameter of hard sphere, a; is a unit vector along 
the line passing through the centers of the spheres at their interaction, B(v — w, uS) = (v — w)u 
is the collision kernel and are the collision factors. F^ are usually assumed to be two 
functionals of /, more precisely speaking, they depend on the density p(t,x) = J R3 f(t,x,v)dv 
at the time t and the point x. 
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In equations (|1.2p and (j 1 . 3 j) . (v,w) and (v',w') are velocities before and after the collision, 
respectively. As for the Boltzmann equation, the conservation of both kinetic momentum and 
energy of two colliding particles gives 

v + w = v' + w', v 2 + w 2 = v' 2 + w' 2 , (1.4) 

which leads to their velocity relations 

v = v — [(v — w)u]u, w = w + [(v — w)u>}uj, (1-5) 

where u G S+. We denote u = v — w, uu = (tiw)w and u± = u — un (see [9] or [H]), thus getting 
another expression of fjl .5[) as follows: 

v = v — u\\, w' = v — u±. (1-6) 

Then the gain and loss terms (|1.2|) and (|1.3|) can be rechanged as 

Q + (f)(t, x, v ) = / F + (f)f(t,x,v — u\\)f(t,x — au>,v — u±)B(u,uj)dujdu, (1-7) 

Q~(f)(t, x, v) = / F~(f)f(t, x,v)f(t,x + aco,v — u)B(u, u))dwdu. (1-8) 

If the factors i ?± are set to be the same positive constant and the diameter a equal to zero 
in the density variables, then the Enskog equation becomes the Boltzmann one that provides a 
successful description for dilute gases. The Boltzmann equation is no longer valid for moderately 
or highly dense gases. As a modification of the Boltzmann equation, the Enskog equation 
proposed by Enskog |10j in 1922 is usually used to explain the dynamical behavior of the density 
profile of moderately or highly dense gases. 

As we know, there are global solutions for the Boltzmann equation in the absence of external 
forces not only in an infinite vacuum but also with large initial data. Existence of such vacuum 
solutions was first considered by Illner & Shinbrot [12] and later by Bellomo & Toscani [3]. 
The global existence of solution was shown by DiPerna & Lions [S] for the Boltzmann equation 
with the large data. There are also some similar results about the Enskog equation without 
any effect of external forces. For example, Polewczak [15] and Arkeryd [1] gave their different 
existence proofs of global-in-time solutions to the Enskog equation without external forces for 
both near-vacuum and large data, respectively. It is worth mentioning that some early works 
on the existence of global solutions to the Enskog equation were given by Cercignani and/or 
Arkeryd (e.g. [2], [3], [5], [6]). Recently, Duan et al. [9] proved the existence and uniqueness of 
a global mild solutions for the Boltzmann equation with external forces in an infinite vacuum 
and many relevant works can be found in the reference. Now there is not yet such similar result 
for the Enskog equation in the presence of external forces. The aim of this paper is to extend 
the result to the case of the Enskog equation, that is, to show the existence and uniqueness of 
such vacuum solution to the Enskog equation (jl.ip with (II. 2|) and (jl.3jl . In Section [2] hypotheses 
of the external forces are given and a Banach space and its operators are constructed. These 
hypotheses are very related to collision of particles of moderately or highly dense gases along 
the bicharacteristic curves and they are much more complicated than those given by Duan 
et al. mentioned above for the Boltzmann equation since the collision operator of the Enskog 
equation is fairly more complicated than that of the Boltzmann equation. In spite of this, the 
two examples shown by Duan et al. in [S] satisfy our hypotheses. Then the estimation of the 
so-called gain and loss integrals is made in Section [3l An existence and uniqueness theorem of 
global mild solution to the Enskog equation in an infinite vacuum is given in Section 0J and some 
remarks on the assumption of the factors F^ are finally made in Section 
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2. Hypotheses and Operators 

In this section we first give some constructive hypotheses of the external forces with the help of 
the bicharacteristic equations of the Enskog equation with these forces and then build a Banach 
space and its operators relative to the Enskog integral equation. 

Let us begin with considering the bicharacteristic equations of the Enskog equation (11,11) 

^§ = V, ^=E(s,X), (X,V)\ s=t = (x,v). (2.1) 

Suppose that such a vector-valued force function E(t, x) allows the above system (|2.ip to have 
a global-in-time smooth solution denoted by 

[X(s;t,x,v),V(s;t,x,v)} (2.2) 

for any fixed (t, x, v) G R + x R 3 x R 3 , and that there exist three functions a«(s; t, x, v) (i = 1, 2, 3) 
such that the solution (|2.2p satisfies the following conditions: 

ax(s;t,x,v) > as s > 0, a\(0;t, x,v) > 0, ct2(s;t, x,v) > 0, as(s;t, x,v) > 0, (2-3) 

X(0; s, X(s; t, x, v) + V(s; t, x, v) — rj) = X(0; t, x, v) + a\(s; t, x, v)r] + q 2 (s; t, x, v)£, (2.4) 

V(0; s, X(s; t, x, v) + £, V(s; t, x, v) — rj) = V(0; t, x, v) - a 2 (s; t, x, v)r] - a 3 (s; t, x, v)£, (2.5) 

either as(s;t,x,v) = or max{ai(s;t,x,v),ct3(s;t, x, v)}/a2(s;t, x,v) < To, (2-6) 

min{(a2(s; t, x, v )) 2 a(s; t, x, v ), a2{s; t, x, v)a(s; t, x, v ), «2 (s; t, x, v)} > ao > 0, (2-7) 

for any s G R+ and (£,r]) G R 3 X R 3 when any point (t,x,v) is fixed in R + x R 3 x R 3 , where 
ao, eo and tq are three fixed positive constants independent of s and (t, x,v), and a(s;t,x,v) is 
denoted by 

a(s; t, x, v) = a'i{s] t, x, v)ct2(s; t, x, v) — a%(s; t, x, v)a' 2 (s; t, x, v), (2-8) 

here a'As;t,x, v) (i = 1,2) represent the derivative with respect to s. For their understanding 
of their hypotheses of the external force in the Boltzmann equation, Yuan et al. [9] took the 
following two examples: E(t,x) = Eo(t) and E(t,x) = e^x + Eo(t) with eo being a positive 
constant. Since the Enskog equation is much more complicated than the Boltzmann one, our 
hypotheses of the external forces for the Enskog equation are much more complicated than those 
given by Yuan et al. mentioned above for the Boltzmann equation. In spite of this, the above 
two examples satisfy the above hypotheses of the external force. These examples are obviously 
suitable for our explanation of the corresponding constructive conditions whether the Boltzmann 
equation or the Enskog one is considered. 

We give the five conditions (I2,3p - (|2.7|) in order to get the following inequalities 

\X(0; s, X(s; t, x,v),V(s;t,x,v) — u\\)\ 2 + \X(0; s, X(s; t, x, v) — au, V(s; t, x, v) — u±)\ 2 

> \X(0; t, x, v)\ 2 + \X(0; t, x, v) + a\(s;t,x, v)u — 002(5; t, x, v)lo\ 2 

and 

| V(0; s, X(s; t, x,v),V(s; t, x, v) — u\\)\ 2 + | V(0; s, X(s; t, x, v) — aw, V(s; t,x,v) — u±)\ 2 

> | V(0; t, x, v)\ 2 + | V(0; t, x, v) — 02(5; t, x, v)u + aas(s; t, x, v)lv\ 2 

along their bicharacteristic curves after collision of particles of moderately or highly dense gases. 
The above two inequalities can be used to control the gain and loss integral terms of the Enskog 
integral equation along the bicharacteristic curves. Therefore this form of the external forces is 
pertinent to collision of particles of moderately or highly dense gases along the bicharacteristic 
curves so that the gain and loss integral terms can be estimated. 

We below give a representation of mild solution to the Enskog equation. Let us first introduce 
a notation /# defined as 

f*(s;t,x,v) = f{s,X(s;t,x,v),V{s;t,x,v)) 



-1 



ZHENGLU JIANG 



for any measurable function / on R + xR 3 xR 3 . Obviously, it can be found that f*(t;t,x,v) = 
fit, x, v) and that / # (0; t, x, v) = f(0, X(Q; t, x, v), V(Q; t, x, v)). 

Along the bicharacteristic curves, the Enskog equation (jl.ip can be also written as 

-^f#(s;t,x,v) = Q(f)*(s;t,x,v), 
as 

which leads to the following integral equation 

f(t, x, v) = / (X(0; t, x, v), V{0; t, x, v)) + I Q(ff{s; t, x, v)ds, (2.9) 

J o 

where fo(x,v) = f(0,x,v). A function f(t,x,v) is called global mild solution to the Enskog 
equation (jl.ip if f(t,x,v) satisfies the above integral equation (|2,9p for almost every (t, x,v) G 
R+ x R 3 x R 3 . 

Then we construct a subset M of a Banach space C(R+ x R 3 x R 3 ), which has the property 
that every element / = f(s, x,v) £ M if and only if there exists a positive constant c such that 
/ satisfies 

\f{t, x,v)\< ch(X(0; t, x, v))m(V{0; t, x, v)) 

where 

h(x) = e- p|x|2 , m{v) = e~ qv \ (2.10) 

for any fixed p and q in (0, +oo). It follows that M is a Banach space when it has a norm of the 
following form 

Il/H = suvilf^x^^h-^X&t^^^m^iViO^x^))}. 

t,x,v 

In particular, 1 1/(0, x, v)\ \ = sup{\f(0,x,v)\h~ 1 (x)m^ 1 (v)}. The initial data /o = f(0,x,v) is 

x,v 

bounded in L X (R 3 x R 3 ). This implies that the total mass of the system is finite. Hence the 
mean free path is sufficiently large if the finite total mass is sufficiently small. This is exactly 
the requirement on the Enskog equation with external forces in an infinite vacuum, which is 
similar to one considered by Illner & Shinbrot [12] for the Boltzmann equation. It is worth 
mentioning that in other cases there are many different classes of functions which can be taken 
as the choice of h(x) and m(v) (see [9], [ll]). For example, in the case of the external forces 
depending only on the time t, one can also choose h(x) = (1 + |x| 2 )~ p and m(v) = e~ qv for any 
fixed p £ (1/2, +oo) and q £ (0, +oo); when p > 3/2, the initial data /o = f(0,x,v) is bounded 
in L X (R 3 x R 3 ); when 1/2 < p < 3/2, the initial total mass might be infinite; this choice of both 
h(x) and m(v) is not suitable for using our method considered in this paper to deal with this 
existence problem of vacuum solutions to the Enskog equation with external forces depending on 
the time and space variables, however, this choice can be applied to the case of the Boltzmann 
equation in the presence of external forces [9]. 

To give global existence, it is necessary to study the properties of the collision operator in 
a Banach space. To do this, by (|1.7|> and (II. 8p . Q(f)#(s;t, x,v) can be first rewritten as the 
difference between the gain and loss terms of other two forms 

Q + (f)#(s;t,x,v)= [ F+(f)f(s,X(s-t,x,v),V(s;t,x,v)-u {] ) 

x/(s, X(s; t, x, v) — aid, V(s; t,x,v) — u±)B(u, u>)du>du, (2-H) 

Q-(f)#(s;t,x,v)= [ F-(f)f(s,X(s;t,x,v),V(s;t,x,v)) 

x/(s, X(s; t, x, v) + auj, V(s; i, x, v) — u)B(u, u>)du>du. (2-12) 



GLOBAL SOLUTION TO ENSKOG EQUATION 5 

Estimation of the collision integrals can be then made by use of a similar argument to that 
developed in the previous work (see [9], [15], [IB]). According to (|2.1ip and (|2.12p . we in fact 
have to estimate the following two integrals: 

I g = / h(X(0; s, X(s; t, x, v), V(s; t, x, v ) — um)) 
Jo JU 3 xSl 

x m(V(0; s, X(s; t, x, v), V(s; t, x, v) — u^))h(X(0; s, X(s; i, x, v) — atu, V(s; t, x, v) — u±)) 

x m(V (0; s, X(s; i, x, v) — auj, V(s; t, x, v) — u±))B(u, uj)dujduds, (2-13) 

Il= f l h{X(0;s,X{s;t,x,v),V(s;t,x,v))) 
JO JR :i xSl 

x m(V(0; s, X(s; t, x, v), V(s; t, x, v)))h(X(0; s, X(s; t, x, v) + au>, V(s; t, x, v) — u)) 

x m(V(0; s, X(s; t, x, v) + am, V(s; t, x, v) — u))B{u, u>)du>duds. (2-14) 

Here I g and Ii are called the gain and loss integrals respectively. Once the estimation of the 
integrals is finished, the global existence result may be obtained by constructing a contractive 
map from a Banach space to itself. Therefore it is one of the best important to estimate the 
above two integrals, especially the gain one. It will be discussed in the next section. 
We finally denote an operator J on M by 

J(f) = f (X(0; t, x, v), V(0; t, x, v)) + / Q(f)#(s; t, x, v)ds. (2.15) 



It will be proved in Section H] that J is indeed a contractive map on a Banach space. This is 
what we need. 

3. Estimation of the Gain and Loss Integrals 

In this section the so-called gain and loss integrals are estimated by use of a similar device to 
one given in [15J for the Enskog equation in the absence of external forces. To do this, we first 
introduce the preliminary lemmas which will be used below. 

Lemma 3.1. Let any z E R 3 , s E R + , (u\\,u±) E R 3 x R 3 with u\\u± = and u> E S 2 with 
uulo > 0. Then 

\z ± sii|| | 2 + \z ± su± =F au>\ 2 > \z\ 2 + \ z ± s(u\\ + u±) =F o-uo\ 2 (3-1) 
for any fixed real number a E R+ . 

Lemma 3.2. Let p > and (z, u) E R x R with Then 



+oo nz 

Vp\u\ 

Lemma 3.3. Let q > 0, -2 < 7 < 1 and z E R 3 . Then 

ur^l-l^<4- + ^. (3.3) 
R3 7 + 2 q i ' 2 

Lemma 3.4. Let any (s, z, u) E R x R 3 x R 3 and u>± be a unit vector perpendicular to u> E S 2 . 
Then 

\z + su + hu\ > \zuj±_ + suu±\ (3-4) 

for any h E R. 

The proof of Lemmas 13.11 and 13.41 is easily given. Lemmas 13.21 and 13.31 can be obtained from 
the transformation of integral variables. 
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Lemma 3.5. Assume that three functions cti(s;t,x,v) (i = 1,2,3) satisfy the external force 
conditions \2. 3\) and {2. 7|), and that 



I l {z 1 ,Z 2 ,t,X,v)= [I | uw | e -p|*i-H»i(«it,x,*)«+aa 3 («;t,x I »)a;| 
JO Jn^xSi 



2 



xe ^q\z 2 - a2 (s;t,x,v)u-aa 3 (s;t,x,v)^ dudu;ds 

for any fixed real number a G R+ . Then 

h{zi,z 2 ,t,x,v) < ii pq (3.5) 

for any [z\, z 2 ) G R 3 x R 3 , (t, x, v) G R+ x R 3 x R 3 , where Ii pq is a positive constant depending 
only on p and q, p > 0, q > 0. 

Proof First let us fix {z\,z 2 ) G R 3 x R 3 and (t, x, v) G R+ x R 3 x R 3 . Put u = a 2 (s; t, x, v)u + 
aa.%(s; t, x, v)u. Then ulo = a 2 (s;t,x,v)(uu)u + aas(s;t, x,v) for uj G S\. By (I2.3p . it thus 
follows that 

Il(zi,z 2 ,t,x,v) < / / \uuj\(a 2 (s;t,x,v))~ A 
Jo Jr 3 xsI 

I t a i (s;t,x,v) _ , aa(s;t,x,v) \2 _,o 

Xe ~P^+ <* 2 (s;t,*,v) u+ a 2 ?. i t,«,«> w l e -fll«-a| dudujds, (3.6) 
where a{s; t, x, v) = a 2 {s; t, x, v) — a\{s; t, x, v)ct3(s; t, x, v). Take r = . Then, by (12. 7p . 

% = (^{ift^x'v))' 2 > ^' ^ replacing the integral variable s with the new variable t and using 
Lemma 13.41 the estimation of the integral on the right side of (|3.6|) thus gives 

r+oo r 

I l (z 1 ,z 2 ,t,x,v) < / / \uuj\(a 2 (s;t,x,v))~ 2 (a{s;t,x,v)y 1 
JO JR 3 xS% 

x ^z^^tuuj^ e -q\z 2 -u\2 dudujdT 

< — f \uuj\ { e- p \ Zl ^ +TU ^ 2 dT\e- q ^-^ 2 ducLu, (3.7) 
«0 jR3 X 5i Uo J 



where wj_ is a unit vector perpendicular to uj. The last inequality in (13. 7p comes from (12. 7D and 
Fubini's theorem. By estimation of the integral on the right side of the last inequality, it then 
follows that 

I l ( Zl ,z 2 ,t,x,i>)<—^= -J_L e -9l*a-«l dudu < + , (3.8) 

where Lemmas 13.21 and 13.31 are used. The proof of Lemma 13.51 is hence completed. □ 

Lemma 3.6. Assume that three functions ai(s;t,x,v) (i = 1,2,3,) satisfy the external force 
conditions \2. 3\) . 112. 6\) and {2. 7\ ), and put 

Ig( Zl ,Z 2 ,t,X,v)= !! \ uu) \ e -p\zi+ ai {s-,t,x,v)u-a a 2(s,t,x,v)^ 
Jo Jn s xSl 

xe -q\z2-a 2 (s^x,v)u+aa 3 (s;t,x,v)u J ^ dudud ^ 

where a > 0. Then 

I g (zi,z 2 ,t,x,v) < I gpg (3.9) 

for any {z\, z 2 ) £ R 3 x R 3 , (t, x, v) G R + x R 3 x R 3 , where I gpq is a positive constant depending 
only on p and q, p > 0, q > 0. 
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Proof. Let us fix (z\, z%) G R 3 x R 3 and (t, x, v) G R+ x R 3 x R 3 . Note that a2(s; t, x, v) > in 
(|2.3p . Put u = un + where uii = a2(s; t, u)wn — aa^{s; t, x, v)lo and u± = a2(s; i, cc, v)u±. 
Then no; = «2(s; t, x, v){uuj) — aa%(s; t, x, v) for uj G Thus 

I g (zi,z 2 ,t,x,v) < / / (|nw| +aa 3 (s;t,x,w))(a2(s;t,2;,t')) _4 

x e _p|Zl Q 2 ( s ;*^^) u_ Q2( s ;*^^) aJ| e _| 3|22-n| dudujds, 

where a(s;t,x,v) = a\{s;t,x,v)) — a\{s;t, x,v)a^{s;t, x,v). To prove this lemma, it suffices to 
consider the case of max{ai(s; t, x, v ), as(s; t, x, u)}/a 2 (s; t, x, v ) < tq in (|2.3p . By repeating a 
similar integral estimation to one given in Lemma |3.5|, the estimate of the integrals on the right 
side of the above inequality gives 

(a) rr r 

I g (zi,Z2,t,x,v) < / / (\uuj\ + aa 3 (s;t,x,v))(a2(s;t,x,v)y 2 (a(s;t,x,v)y 1 

JO JR^xS* 

x e - p ^^ +T ^^ 2 e- q \ Z2 -^ 2 dudujdT 

(b) I/ 1 ( r+oo -\ 

< — / l\uco\ / e- p \ z ^^ +T ^^dT + aT%\e- q \ z *-^dudw 
«o Jn?xs% I Jo J 

^ 1 f ( \ulo\ 2 ) fi|3 , , ( d ) 4vr 2 /^F 2 \ / 4 1 x 

< — / + «r 2 e -«l«»-«l < V_ + ar 2 - + — , 3.10 

for a3(s; t, x, u) 7^ 0, where lo± is a unit vector perpendicular to uj, (a) is obtained by first making 
the transformation r = and then using (|2.6I) and Lemma 13.41 (b) is given by (12.70 . (c) 

is obtained by Lemma 13.21 and (d) results from Lemma 13.31 This hence completes the proof of 
Lemma [3^61 □ 



By Lemmas 13.51 and 13.61 we can hence give the estimates of the gain and loss integrals as 
follows. 

Lemma 3.7. Let I g and I\ be the same integrals as defined by \2.13\) and \2.1J$ , respectively. 
In the integrals, a is a positive constant. Suppose that all the five external force conditions 
112. 3\) - (277\ ) hold for these functions X(s; t, x, v) and V(s; t, x, v) defined by the solution 12. 2\) to 
the system \2. and that h(x) and m(v) are the same as in 12.10\) . Then it follows that 

I g < Kh(X(0; t, x, v))m(V(0; t, x, v)), (3.11) 

// < Kh(X(0; t, x, v))m(V(0; t, x, v)), (3.12) 
for any (t,x,v) G R+ x R 3 x R 3 and some positive constant K. 

Proof. Let us first estimate the loss integral. By using (|2.4p and (|2.5|) . the loss integral (|2,14p 
can be rewritten as 

J l — f f \ uu \ e ~p\X(0;t,x,v)\ 2 e -p\X(0;t,x,v)+ai(s;t,x,v)u+aa2(s;t,x,v)uj\ 2 

JO iR 3 xS| 

x e -ql^( ; ^ . I ^' u )| 2 e - | ^|V'(0;t,x,^))-a2(s;^,a;,^))u-aa3(^i;^,x,■!;)a;| 2 ^ ^^ ^ w ^ s 

= h(X(0; t, x, v))m(V(0; t, x, w))J/(X(0; t, x, v), V(0; t, x, v),t, x, v)). (3.13) 
It follows from Lemma 13.51 that 

/, < h{X(0;t,x,v))m(V(0;t,x,v))I lpq . (3.14) 
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Then the estimation of the gain integral will be made below. Similarly, by using (|2.4p and 
(|2.5p . the gain integral (|2.13p is as follows: 

t f f I \ —p\X (0;t ,x ,v)+a\ (s;t,x,v)uu \ 2 —p\X(0;t,x,v)+ai (s;t,x,v)uj_—act2 (s;t,x,v)ui\ 2 

/o J-R 3 xSi 



x e -<l\V(®\t,x,v)-U2{s;t,x,v)u\\ 2 e ~q\V{0;t,x,v)-a 2 {s;t,x,v)u A _+aa 3 (s;t,x,v)w\ 2 (3.15) 

Using Lemma 13. 11 we have 

j < e -p\X(0;t,x,v)\ 2 e -q\V(0;t,x,v)\ 2 f f ^ u ^ e - P \X(0;t,x,v)+a 1 (s;t,x,v)u-aa2(s;t,x,v)u)\ 2 
9 ~ JO JRSxSl 

x e — <?l V 7 (0;t,x,v)-a 2 (s;t,x,v)u+aa 3 (s;t,x,v)w\ 2 dudujds 

= h(X(0; t, x, v))m(V{0; t, x, v))I g (X(0; t, x, v), V(0; t, x, v),t, x, v). (3.16) 
It follows from Lemma 13.61 that 

I g < h(X(0; t, x, v))m(V(0; t, x, v))I gpq . (3.17) 
The proof of Lemma 13.71 is hence completed. □ 

4. Existence and Uniqueness 

In this section we show a result about the existence and uniqueness of such vacuum solution 
to the Enskog equation in presence of external forces. To do this, we first define a set Mr by 

M R = {f: Il/H <RJeC(H + xR 3 xR 3 } (4.1) 

and then assume that i ?± are two functionals on Mr such that the so-called locally Lipschitz 
condition 

\F±{f)-F±{g)\<L{R)\\f-g\\ (4.2) 
holds for any f,g € Mr where Mr is defined by (|4.ip and L{R) is a positive nondecreasing 
function on R + . Thus we can get the following lemma. 

Lemma 4.1. Suppose that all the five conditions \2. 3\) -i2.7) hold for any X(s;t,x,v) and 
V(s; t, x, v) defined by the solution \2. 2\) to the system \2.1\) . and that the factors in the 
collision integrals Q^{f)^{s;t,x,v) defined by \2.1V\) and \2. 12\) are two functionals satisfying 
the inequality Let h(x) and m(v) be the same as in \2. 10\) . In the collision integrals, a is 

a positive constant. Then the following inequalities hold: 

' \Q + (f)*(s;t,x,v)\ds < C(R)h(X(0;t,x,v))m(V(0;t,x,v))\\f\\ 2 , 

1 \Q-(f)*(s; t, x, v)\ds < C(R)h(X(0; t, x, v))m(V(0; t, x, v))\\f\\ 2 
for any f £ Mr, where C(R) is a positive nondecreasing function on R + . 

Proof. It can be first found from the assumption (14. 2p of the two functionals that there 
exists a positive constant L(R) = L{R)R + |F+(0)| + |F~(0)| such that ^(f)] < L(R) for any 
/ £ Mr. It follows from (12TTT1) and ([2/T2"]) that 

-{f)*{s;t,x,v)ds<L{R) f [ \\f\\ 2 h(X(0;s,X(s;t,x,v),V(s;t,x,v)-u ll )) 

Jo Jr 3 xs| 

x m(V(0; s, X(s; t, x, v), V(s; t, x, v) — u\\))h(X(0; s, X(s; t, x, v) — au, V(s; t, x, v) — u±)) 

x m(y(0; s, X(s; t, x, v) — au, V(s; t, x, v) — u±))uudududs, (4-3) 

Q-(f)#(s;t,x,v)ds < L(R) f [ \\f\\ 2 h(X(0;s,X(s;t,x,v),V(s;t,x,v))) 

Jo Jr s xsI 

x m(V(0; s, X(s; t, x, v), V(s; t, x, v)))h(X(0; s, X(s; t, x, v) + au, V(s; t, x, v) — u)) 
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x m(V(0; s, X(s; t, x, v) + aw, V(s; t, x, v) — u))uu>duiduds. (4-4) 
By (ffTTT) and (f3~T2l . (g^D and (JO]) give 
r-t 

Q + (/) # (r, x, w)dr < L(R)Kh(X(0; t, x, v))m(V(0; t, x, v) 

t 

Q~(f)*(T,x,v)dT < L(R)Kh(X(0;t,x,v))m(V(0;t,x,v 

Take C(R) = L{R)K. It obviously follows that Lemma 14 . 1 1 holds . □ 
Then we can get the following theorem. 



Theorem 4.2. Suppose that all the five external force conditions $2.3\) -i 2.7) hold for these 
functions X(s;t,x,v) and V(s;t,x,v) defined by the solution \2.2\) to the system 12. and 
that the factors in the collision integrals Q ± (f)(t,x,v) defined by M.2\) and are two 

functionals satisfying the inequality I n the collision integrals, a is a positive constant. 

Then there exists a positive constant Rq such that the Enskog equation M.l\) with M.2\) and 
hl. 6 J\) has a unique non-negative global mild solution f = f(t, x, v) E M# through a non-negative 
initial data fo = fo(x,v) when 

sup{/ (X(0; t, x, v), V(0; t, x, v))^ 1 (X(0; t, x, v)) m - 1 (V(0; t, x, v))} 

t,x,v 

is sufficiently small, where h(x) and m(v) are the same as in \2.1(K) . 

Theorem 14.21 shows that there exists a unique global mild solution to the Enskog equation 
(jl.ip given by (|l-2j) and (jl .3f) with the initial data near vacuum if a suitable assumption of the 
external force is given. As in p~3), we below give our proof of Theorem 14.21 

Proof. By (|2.15p and Lemma 14.11 we have 

| J^h- 1 (X(0; t, x, v))m~ l (V(0; t, x, v)) 
< |/o(X(0; t, x, v), V(0; t, x, v^h-^XiO; t, x, v)) m - 1 (V(0; t, x, v)) + 2C(R)\\f\\ 2 

< R/2 + 2C{R)R 2 

for any / G Mr and fo with ||/o|| < R/2. Since C{R) is a positive nondecreasing function on 
R + , it follows that || J(/)|| < R for sufficiently small R. Therefore J is an operator from Mr to 
itself for sufficiently small R. Similarly, it can be also found that J is a contractive operator on 
Mr for some suitably small R. Thus there exists a unique element / G Mr such that / = J(f), 
i.e., (12. 9p holds. It then follows from the same argument as the one in [12] (or see [H], [T9j ) that 
if fo(x,v) > then f(t,x,v) > 0. Hence the proof of Theorem 14.21 is finished. □ 

5. Remarks on the Assumption of the Factors F^ 1 

In this section we make some remarks on the locally Lipschitz assumption (14. 2D of the factors 
of i ?± appearing in Theorem 14.21 given in the previous section. 

We begin with a different kind of locally Lipschitz condition of F . It was originally given by 
Polewczak [15 1 as follows: 



\F ± (f)-F ± (g)\<L (R) 



/ f(t,x,v)dv - / g(t, 

JR3 7R3 



x, v)dv 



(5.1) 



holds for any / = f(t,x,v),g = g(t,x,v) £ Mr where Mr is defined by (|4.ip and Lq(R) is a 
positive nondecreasing function on R + . Note that assumptions (|2,4p and (|2.5p have the following 
properties: 

X(Q; t,x + £,v)- X(0; t, x, v) = a 2 (s; t, x, v)£ (5.2) 

and 

V(0; t, x, v + rj) — V(0; t, x, v) = «2(s; t, x, v)r] (5-3) 
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for any (£, 77) G R 3 x R 3 when any point (t, x, v) is fixed in R + x R 3 x R 3 . By (|5.2p and (|5.3 
we have 

<9X(0; t, x, v) dV(0;t,x,v) 



thus giving 



dx dv 
dX{0;t,x,v) dV(0;t,x,v) 



a 2 (s;t,x,v), (5.4) 
> a > (5.5) 



9a; 9f 

because of assumption (|2.7j) . Put L(-R) = Lq(R) J r3 m(v)dv / a^. Then, by use of (15. 5p . (15. 1|) 
gives (|4.2p . This means that (|5.ip is a stronger assumption than (|4.2p when the external forces 
satisfy the assumptions (|2.3p - (|2.7p . Therefore the locally Lipschitz assumption (|4.2p is at least 
mathematically very useful in a more general case. 

Now let us recall the Enskog equation for our further understanding the locally Lipschitz 
assumption (|4.2p of the factor F under the assumptions (|2.3p - (|2.7p of the external forces. The 
Enskog equation can be roughly divided into two classes: the standard and the revised one. In the 
standard Enskog equation (see [7j, [10j . |15]). i ?± are defined by a geometrical factor Y which 
is a contact-value pair correlation function of the hard-sphere system at uniform equilibrium 
and depends on the density p(t,x), i.e., they are given by F + = Y(p(t,x — aw/2)) and F~ = 
Y(p(t, x + atv/2)). For a fairly rare uniform gas of one particle with mass m, it can be found in [7] 
that the value of Y is approximatively expressed by Y(p(t, x)) = [1 — llbp(t, x)/8]/[l — 2bp(t, x)] 
where b = 27ra 3 /(3m). It can be easily shown that in this case the factor F^ = Y satisfies (15. ip 
and is locally Lipschitz as defined in (14. 2p with the external forces of the assumptions (I2.3p - (|2.7|) 
for R G (0,Rq], where Rq is some suitably small positive constant. Generally, the dependence 
of the function Y on the local density p(t, x) is of the form 

+00 

Y(p(t, x)) = l + Y J k[2Tra?p{t, x)/3]\ (5.6) 

i=l 

where b% are given in terms of the virial coefficients Bi appearing in the equation of state for the 
hard sphere system. We cannot know whether the series (15.61) converges when one of the two 
different locally Lipschitz conditions (14. 2p and (15. ip is satisfied. Even if this series converges, 
we cannot yet know whether one of the assumptions (|4.2p and (I5.ip of F holds for any factor 
Y of the above form. Of course, if F ± = Y is a factor defined by the form (15. 6|) with finite 
terms, then (15. ip holds and so F satisfy (I4.2|) when the external forces are of the assumptions 

In the case of the revised Enskog equation (see pQ , [16] , [20] ) , F ± are expressed by a contact- 
value pair correlation function G of the hard-sphere system at non-uniform equilibrium. The 
form of G is given by the Mayer cluster expansion in terms of local density p(t, x). The function G 
depends on the position x, the vector x±auj and the density p(t, x), i.e., F + = G(x, x—auj, p(t, x)) 
and F~ = G(x, x + au, p(t, x)). In fact, G(x, y, p(t, x)) = exp(— (3$(\x — y\))G(p(t, x)), where (3 
is a positive constant, <£(|x — y\) is a potential of two interaction spheres at the positions x and 
y, and G is a functional of the following form [16] 

G(p(t,x)) = 1 + J F(12|3)p(3)cix3 + i J J K(12|34)p(3)p(4)di3(fe4 



//■"/ F(12 ' 34 " " " k )P^P^ • • • P( k ) d ^dx^ ■ ■ ■ dx k + • • • , (5.7) 

here p{k) = p(t,Xf~) and V(12|34- • • k) is the sum of all the graphs of k labeled points which 
are biconnected when the Mayer factor exp(— /33>(|xj — Xj\)) — 1 are added. In contrast to 
the standard Enskog equation, the revised Enskog equation possesses an H-function [T7]. It 
can be also known that the revised Enskog theory for mixtures is consistent with irreversible 
thermodynamics, including Onsager reciprocity relations (see [21] , [22], |23|). But we cannot 
know whether the above series (|5.7p converges under the assumption (|4.2p or (|5.ip . We cannot 
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yet know whether one of the assumptions (14. 2p and (|5.ip of F holds for any functional G of 
the above form. 

The assumption (14. 2p or (15. ip of is satisfied by some geometric factors present in the 
truncated Enskog equations in both standard and revised cases. In the standard case the 
geometric factor Y considered above can be assumed to be of the truncated form (I5.6|) with finite 
terms while in the revised case an obvious example is that the functional G can be truncated 
to be a positive constant. It can be found that these factors satisfy the two assumptions of 
F^ when one assumes that both L(R) and Lq(R) are two positive constant functions on R + . 
Therefore the above two assumptions are completely suitable for our understanding evolutions 
of moderately or highly dense gases by use of our investigation of the properties of the Enskog 
equation. 

Acknowledgement. The author would like to thank the referees of this paper for their valuable com- 
ments and suggestions on this work. 
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